Odd-frequency Cooper pair amplitude around a vortex core in a chiral p-wave 

superconductor in the quantum limit 
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Solving the Bogoliubov-de Gennes equation, we study the spatial structure of odd-frequency s- 
wave Cooper pair amplitudes at each quantized energy level of vortex bound states in chiral p-wave 
superconductors. For zero energy Majorana states, the odd-frequency s-wave pair amplitude has the 
same spatial structure as that of the local density of states even in atomic length scale. This relation 
also holds in finite energy bound states for single vortex winding anti-parallel to the chirality, but 
not for parallel vortex winding. The double winding vortex case is also studied. 
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I. INTRODUCTION 

Studies of spatial-variation in unconventional super- 
conductors give valuable information to know characters 
of unconventional superconductivity. Among these stud- 
ies, such as for ferromagnct/superconductor junctions,— ~— 
normal metal/superconductor interfaces^—, and vortex 
states ; 10 ' 11 it has been shown that their low energy quasi- 
particle states are tightly related to the induced odd- 
frequency Cooper pair amplitude.— Nowadays, it is ac- 
cepted that the Andreev bound stated and anomalous 
proximity effect^ are interpreted as the emergence of 
the odd-frequency Cooper pairs.— 

Since quasiparticles in superconducting states have to 
satisfy the anti-commutation relation, the pair potential 
must change its sign after a permutation of two quasi- 
particles is done. The minus sign comes from the combi- 
nation of orbital components (odd or even parities), spin 
component (singlet or triplet) and Matsubara frequency 
(even or odd frequency). Conventional even-parity spin- 
singlet pairing and odd-parity spin-triplet pairing are 
classified to the even-frequency pairing. On the other 
hand, within the odd-frequency pairing, even-parity spin- 
triplet pairing^ and odd-parity spin-singlet pairing^ are 
also possible. Here, we study odd-frequency pair ampli- 
tudes induced around a vortex in even-frequency pair- 
ing supcrconductorsi 10 ' 12 When spin-flip mechanism is 
absent for quasiparticles, an odd-frequency spin-singlet 
odd-parity (OSO) pair amplitude appears around vortex 
of even-frequency spin-singlet even-parity (ESE) super- 
conductors, or an odd-frequency spin-triplet even-parity 
(OTE) pair amplitude appears around vortex of even- 
frequency spin-triplet odd-parity (ETO) superconduc- 
tors. We discuss the latter case in this worki^S 

Most of previous studies for odd-frequency pairing 
were based on quasiclassical Eilenberger theory which 
is valid in the limit 2(fcp£)~ 1 = A/Ep <C 1 for super- 
conducting gap A and Fermi energy Ep. kp is Fermi 
wave number and £ is superconducting coherence length. 
When the limit A/Ep 1 is not satisfied in the strong- 
coupling superconductors or we focus fine energy struc- 



tures, low-energy bound states around a vortex core are 
quantized to Caroli-de Gennes-Matricon (CdGM) states 
with energy splitting of the order A 2 /^pi 17 ' 18 The quan- 
tum nature associated with CdGM states has been di- 
rectly observed through the STM-STS in an anisotropic 
superconductor— and the quantum depletion of the par- 
ticle density in Fermi gases^S To study the quantum limit 
case within A ~ Ep , we have to consider the eigen-states 
by solving the Bogoliubov-de Gennes (BdG) equation, 
without using the quasiclassical approximation. 

Recently vortex bound states in chiral p-wave su- 
perconductors attract much attention, because Majo- 
rana state appears at exactly zero energy in the bound 
states.— ~— In a Majorana zero mode of the half quan- 
tum vortex, particle and hole states are equivalent each 
other, giving rise to the non-Abclian statistics of its host 
vortices^ Note that in the case of an integer vortex 
state with Majorana zero modes, the statistics is non- 
trivia h 27 ' 28 On the other hand, the quasiclassical study 
in superfluid 3 He and superconductors showed the rela- 
tion between odd-frequency pair amplitudes and the local 
density of states (LDOS) in Majorana bound state^^ 
Most recently, the issue on the relation between Majo- 
rana zero modes and odd frequency pairing has been 
addressed in a quantum nanowire.— The structures of 
bound states are distinguished by the direction of vortex 
winding, i.e., parallel or anti-parallel to the chirality of 
chiral p-wave superconductivity^ Within the quasiclas- 
sical theory, it has been revealed that the zero energy 
DOS in the parallel (anti-parallel) vortex is fragile (ro- 
bust) against nonmagnetic impurities ; 32 ' 33 which is un- 
derstandable with the odd-frequency pairing localized at 
the vortex^ Hence, since the zero energy state has many 
facets, it is important to study how the odd-frequency 
pair amplitude is related to the bound states, including 
Majorana state, in chiral p-wave superconductors. 

The purpose of this paper is to clarify the relation 
between the CdGM states and the odd-frequency pair- 
ing. To capture the quantum limiting behaviors inside 
the core, we here utilize the BdG equation which is ca- 
pable of describing the rapid oscillations of wave func- 
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tions in atomic length scale of fcp 1 and exactly deals 
with the Majorana mode. The studies in the quantum 
regime have not been done with the quasiclassical the- 
or y 10 ' 11 ' 34 because of the lack of the fcp 1 -scale physics. 
We also note that in previous works, such as Refs.1lH8ll2"9l. 
odd-frequency pairings have been discussed in connection 
with surface bound states. However, the vortex bound 
state which we consider here may exhibit distinct behav- 
iors from surface one in the sense that each level of the 
CdGM states is well isolated in the scale of A 2 /Ep. In 
addition, the quasiparticle structure inside the core is de- 
termined by not only the bulk topology but also the vor- 
tex winding number £2r£L Hence, we here unveil how the 
vortex winding number affects the odd-frequency pairing 
around the core. 



II. FORMULATION OF BOGOLIUBOV-DE 
GENNES EQUATION 

We consider a cylindrical system of isotropic chiral 
p + wave superconductor with radius R = 200fcJ 1 and 
two-dimensional isotropic Fermi surface i 23 ' 24 ' 31 There- 
fore, the cylindrically symmetric order parameter with 
vortex winding w € Z is given by 

A(r,k) = A + (r)e iwt) Y + (k) + A_ (r)c i(t0+2)e Yl (k), (1) 

at the position r = r(cos#,sin(9), where the pairing func- 
tion is given by Y±(k) = (k x ± iky)/\/2k-p for relative 
momentum k = {k x , k y ). In the p + -wave superconductor, 
the Cooper pair has internal angular momentum L z = 1 
along z-axis. In the right-hand side of Eq. ([T]), the first 
term is the dominant ^-component with a vortex, and 
the second term of p_ wave indicates the small induced 
component of opposite chirality around the vortex. 

The cigen-energy Eq and wave functions u q (r), u q (r) of 
quasiparticles are calculated by the BdG equatio n 23 ' 24 ! 31 



ffo(r) 
-H*(r) 



H(r) 
-ffo(r) 



= E a 



u q (r) 

Uq( r ) 



(2) 



where H (r) = (<9 2 + <9 2 )/2A/ - E F with mass M, and 
n ( r ) = -| E m {^ (m) > A m (r)} with in = zb and = 
^fie ±l6 (d r ± ir~ 1 dg)/ \f2k-p in chiral p-wave superconduc- 



tors. k F = (2ME F 
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Throughout this letter, we use 



units S = /cb = 1- For the cylindrical symmetric system, 
wave functions are given by 

Uq (r) = u q (r)e^ e , « q (r) - v^e^-^ . (3) 

The eigen-states are labeled by q = (qe,v), where v as- 
signs eigcn-states of the BdG equation for real functions 
M q (r) and Vq(r). To numerically solve the BdG equation, 
we utilize the Bessel function expansion method, where 
Uq(r) and u q (r) are expanded in terms of the Bessel func- 
tions J q (r) as 



« q M = £c>? s) (r), S q (r) 



3 



where the normalization condition imposed on Uq(r) and 
v q (r) reduces to 27r^(C 2 + D 2 ) = 1. Here, we in- 
troduce the set of the orthogonal functions, tp^ ( 



V2 



^ +l( ra J ^H where k > 



(9) 



and a) is the 



R "J 

j-th zero of J q {r). This reduces the BdG equation ([2} to 
a matrix eigenvalue problem i 18 ' 35 ' 36 

For simplicity we assume that the pairing interaction 
g p works only for chiral p-wave components and that the 
quantization axis of spin is parallel to the d-vcctor of 
spin-triplet pairing. This implies a singular vortex state, 
accompanied by spin-degenerate Majorana zero modest 
Although in the half quantum vortex with a non-Abelian 
Majorana mode the d-vector is transverse to the spin 
quantization axis, the low energy quasiparticles are com- 
monly describable with the BdG equation Here 
we treat Zeeman energy to be negligible at enough low 
fields. Then, the selfconsistent condition A±(r) for chiral 
p-wave pair potentials is given with the imaginary part 
of the retarded Green's function T by 

A±(r)= 3p ^±,tri P ict(r, £)/(£), (5) 

\E\<E mt 

where r denotes the center-of-mass coordinate. Here, 
f(E) denotes the Fermi distribution function and 
even- frequency p±-wave pair amplitudes are given by 

,triplet 

(r, E) = (J±. n (r, E) + F±, it (v, E))/2 with 



F m ^{r,E) 



q 



[u q (ri)u*(r 2 )] 



xS(E - E, 



q 



(6) 



using Vy^ obtained from , p( m ) with r— s-ri2 = r 2 — ri . 
•FmM is obtained from Eq. (|6|) with replacing (u q , w q , Eq) 
to —Eq). It is found that for ETO superconduc- 

tors without a magnetic Zeeman term, even-frequency 
spin-triplet p±-wave pair amplitude J 7 ±^(r, E) satis- 
fies F±^^(v,E) = -J±,tJ,(r, — E) and j±^(r,£?) = 
F± m {t,E)& The BdG equation © with Eqs. © and 
© gives a closed set for the self-consistent calculation of 
A±(r) and wave functions of the eigen-states. 

From the selfconsistent solutions, we calculate the 
LDOS with Uq = Uq(r) and u q = i; q (r) as 

N(r, E)=J2 [SI" S(E - Eq) + \vq\ 2 8{E + Eq)] . (7) 
q 

The odd-frequency spin-triplet s-wave pair amplitude is 
given by F s . triplet = {T s ^i+T s ^)/2 with 



F,, n ( T , E)=J2 Vq(r)uq(r)S(E - Eq) 
q 

7- s4t (r, E)=Y, Uq(r)vq(r)S(E + Eq). (8) 



Note that e l< - w+1 ^J 7 s ,4,t(r, E) can be real, since u q (r) 
and Vq(r) are obtained as a real function from Eq. © 
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FIG. 1: (Color online) (a) Profile of self-consistent chiral 
p-wave pair potential |A±(r)| around a vortex as a function 
of radius r from vortex center. Density plot of (b) N(r, E) 
and (c) -^s, triplet (r, E) with 8 — 0. Left (right) panels are for 
the parallel (anti-parallel) vortex winding. 



with an appropriate U(l) phase of A±. Since 
ETO superconductors without a Zeeman field hold 

F*,u( T > E ) =:F *,u( T >- E ) and Fs,u( T ' E ) = jr s ,lt( T > E )' 
spin-singlet components in ETO superconductors, such 

as Fgrf^TfE) — Fg^frfE), vanish^ Equation (JSJ is 
scaled so that max{e"^ +1)e j;,t r ipiet(r,£ = 0)} = 
max{./V(r, E = 0)}. In figures, the LDOS and pair ampli- 
tude are normalized by these maximum values at E~Q. 



III. NUMERICAL RESULTS 

In our calculation, we set g p = 1.65, E cu t = 2.0EF, T = 
0, so that A~0.12?f when w = 0. First, we consider two 
cases of the parallel (w = 1) and anti-parallel (w = —I) 
vortex state. The obtained pair potentials |A±(r)| are 
shown in Fig. QJa). Around the vortex core of A+(r), 
the opposite chiral component A_(r) is induced. The 
induced component for w = 1 is smaller than that for 
w = — 1, since the winding of A_(r), w + 2 = 3, is larger 
than w + 2 = 1 in the anti-parallel case. 

We examine the relation of LDOS N(r,E) and odd- 
frequency pair amplitude .Fa,tripiet(r, E) near the vortex 
core, which are shown in Figs. |TJb) and QJc) as a func- 
tion of r and E. There, we see the quantization of CdGM 
bound states on E, and the quantum oscillations of wave 
functions with respect to r. For parallel vortex winding, 
the LDOS N(r, E) is not symmetric for E <-> — E, since 
the low energy LDOS at r = is finite only at positive 
energy as shown in the left panel of Fig. (Ub). On the 
other hand, the odd- frequency J r s,tripict(r, E) is symmet- 
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FIG. 2: (Color online) r-dependence of LDOS N(r,E), 
7" a ,triplet(r, E), and J±,tri P iet(r, E) at 6 = 0; (a) Zero-energy 
Majorana bound state, and (b) the second and (c) third lowest 
energy bound states. Left (right) panels are for the parallel 
(anti-parallel) vortex winding. 



ric for E-^—E^ and shows oscillation between positive 
and negative values as a function of r, as shown in the 
left panel of Fig. QJc). It is seen from the right panels 
of Figs.[TJb) andfljc) that for anti-parallel vortex wind- 
ing, .Fs,tripiet(r, E) has the spectral evolution similar to 
N(r,E). Both arc positive and symmetric for E <H> — E 
at \E\<A. 

To examine whether the odd-frequency pair amplitude 
has the spatial structure exactly same as the LDOS or 
not, we compare the profiles as a function of r for each en- 
ergy level. Figure^a) shows N(r, E) and J" Sj 

triplet 

•?p, triplet (r, E) for Majorana bound state at E = 0. 
Here, we confirmed that the relations N(r, E — 0) cx 
^.triplet (r,-E = 0) and 



E ±i triplet (r,E = 0) = 0, 



(9) 



at 9 = in the fully quantum level of fc F , by solving the 
BdG equation (J2J). This implies that the LDOS of the 
Majorana state consists of odd-frequency pairing ampli- 
tudes without including even-frequency ones. Substitut- 
ing the Majorana condition it q (r) = u* (r) into Eqs. ([7]) 
and @ at E~0, we obtain exactly 

N(r, E = 0)oc c^ w+1 ^ e E sMpict (r, E - 0). (10) 

In the zero-energy state for w — 1, N(r,E = 0) = 
-7 r s,tripict(r, E = 0) = at the vortex center r = 0, because 
both itq(r) with qg = 1 and «q(r) with qg — w — 1 — — 1 have 
non-zero winding. For anti-parallel winding w = — l, since 
qg =0 both for u and v, N(r, E = 0) and ^.triplet ( r i E = 
0) have peak at r = 0. The difference by the vortex wind- 
ing direction is due to the angular momentum of CdGM 
vortex bound states. 
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The relation of the LDOS and the odd-frequency pair 
amplitude was studied by the quasiclassical theory for 
surface bound states in superfluid 3 He and superconduc- 
tors However, since the rapid oscillation of the length 
order kp 1 is factorized out by the quasiclassical approxi- 
mation, the theory discusses spatial variations only in the 
length scale of superconducting coherence length. This 
corresponds to the behavior of envelop function of the 
oscillating behavior in Fig. H Note that the STM-STS 
with a O.lnm spatial resolution at low temperature*^ 
has succeeded in unveiling quantum limiting behaviors 
of a vortex core, whose microscopic structures are well 
understandable with the BdG theory but not with the 
quasiclassical theory. 

Next, we study bound states at Ey^O in Figs.[2jb) and 
GJc). In the parallel vortex winding case (left panels), 
the second lowest energy state at E = 0.019Ep has qg = 
for u and qg — w — 1 = -2 for v. Thus, N(r,E) ~ |w| 2 
has a sharp peak at the vortex center r = and the 
odd-frequency pair amplitude .Fa triplet (r, E) ~itv* cx e l2S 
vanishes at r = 0. In Figs. HJb) and[2Jc), J 7 ,. triplet (r, E) 
also shows oscillation between positive and negative val- 
ues at r > 0, implying N(r, E) ^ J" s ,tri P iot(r, E). The third 
lowest energy state at i? = 0.032I?F has qg = — 1 for u and 
qg — w — 1 = —3 for v. Also in this state, we see that 
iV(r, E) ^ J-^, triplet (r, E). For the even-frequency pairing 
amplitude, we see that .F+ , triplet (r, E) ^ 7"- .triplet (r, E). 

For the anti-parallel w = — 1, as shown in the right 
panels of Figs. EJb) and EJc), the second lowest en- 
ergy state at E = 0.0I6I?f has qg = 1 both for u and 
v. The third lowest energy state at E = 0.029.Ef has 
qg = 2 for u and v. Even in these cases, we find that 
N(t,E) oc e-^+^J^tripietO,^ approximately holds, 
that is, the relation in Eq. (TTUJ) can be extended to a finite 
E. This comes from the symmetric LDOS structure for 
E <H- — E and the ^-independent structure w?;*c>ce when 
w = — 1, since the BdG equation ([2]) holds the particle- 
hole symmetry, (u, v) -H> (v* , u*) for E^—E. We also see 
that J r +i tripiet(r,-B)~-7 r -,tri P iot(i", E), but there are small 
deviations between them. 

To study the differences due to vortex winding direc- 
tions in other energy states, in Fig. Eta), we present the 
^-dependence of N(r,E) and e -^ w +^ e J" s , t ri P iet(r, E) at 
r = 2k" 1 . We find er<*>+W .F aitriplot (r, E) £ N(r,E) 
for parallel vortex winding, except for E — 0. Only at 
E = 0, e _t ^ +1 ^J r S) tripiet(r, E) has the same peak struc- 
ture as that of N(r,E), while it has negative value at 
E 0. In the anti-parallel vortex winding case, we 
confirm e - ^ u,+1 ) 9 .Fg ) tripletO r ) E)ocN(r, E) for any bound 
states with \E\ < |A|. For the scattering state at \E\ > |A|, 
however, this relation does not hold. 

The ^-dependences of odd-frequency triplet (r, E) 
and even-frequency T± triplet ( r ) E) arc compared in 
Fig. E£b). From the symmetry relation for E •<->• 
—E, J" s ,triplet (-^i .triplet) is tne even " (odd-) function 
of the real energy E. For anti-parallel vortex wind- 
ing, we see that J"+,tri P let(r, E) ~ .F-,tripiet(r, E). How- 
ever, this is not satisfied for parallel vortex wind- 
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FIG. 3: (Color online) (a) U-dependence of N(r, E) and 
■^.triplet (r, E) at r — 2k^ 1 and 6 — 0. (b) Independence of 
^ r s,tri P ict(r, E) and T± .triplet (r, E). We use Lorenz function of 
half-width 0.001-E'f for 5(E) here. Left (right) panels are for 
the parallel (anti-parallel) vortex winding. 




FIG. 4: (Color online) Density plot of (a) N(r, E) and (b) 
J^, triplet (r, E) for the double winding vortex with w = —2. 
(c) r-dependence of |A+(r)| and |A_(r)|. (d) r-dependence 
of LDOS N(r,E), J" s ,tri P iet(r, E), and j±,tri P iet(r, E) in the 
lowest energy state at E — 0.005^. We set 8 — 0. 



ing. The odd- frequency pair potential A Sj tripiet(r) oc 
/ -^s, triplet (r, E)f(E)dE has to vanish, as confirmed in 
the formulation on Matsubara frequencies. Thus, for 
•^s, triplet (r, E). positive contributions at low energy are 
canceled by the negative contribution at higher energies. 

Lastly, we study the case of a double winding vor- 
tex. Here, zero-energy Majorana states do not ap- 
pear, since the energy levels of the vortex bound 
state for w — ±2 have a small gap of the order 
A 2 /^' 24 ' 36 ' 38 Therefore, this vortex does not exactly 
holds e-^+^J^tripietfr-EO ^ N(r,E) inside the vor- 
tex core. However, in vortex bound states for w = —2, 
N(r,E) in Fig.gfa) and J" S)t ripiet(r, E) in Fig.gfb) have 
similar structure each other. In the chiral p-wave pairing, 
since the induced p- component shares the core region 
without winding (w + 2 = 0), there appears the chiral 
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domain wall between outer p + and inner p_ regions at 
r ~ 19/cp 1 ~ £ as in Fig. g^c)^ Thus, while the low 
energy bound states for w = —2 appear far from the vor- 
tex center, they reduce to the bound states at the do- 
main wall. The quasiparticle bound at the domain wall 
has similar structure to surface Majorana state in chi- 
ral p-wave superconductors^ It is seen from Fig. EJd) 
that the relations e- l{w+1)9 J r sM - ip i ct (r,E)(xN(r,E) and 
-7 r +,tri P ict(r,-B) ~ 7"_,tri P iet(r)-E) approximately hold in 
the low energy, while there are small deviations. This is 
similar behavior to that of w = — 1 in Fig. [2] 

The LDOS N(r,E), which we present in Figs. 1-4, is 
detectable through a STM-STS experiment. Since it is 
demonstrated in Figs. 2 and 3 and Eq. (fTQ|) that the odd- 
frequency pairing amplitude, Js, triplet (r, E), traces the 
spatial shape of the zero energy LDOS, the characteris- 
tic behaviors may be unveiled by using superconducting 
STMfi^ The local Josephson current between the super- 
conductor and the superconducting STM tip is allowed 
only when the symmetry of local pair amplitudes in the 
superconductor is matched with that of the supercon- 
ducting STM tip. Hence, the sharp peak of the odd- 
frequency s-wave pair amplitude around the vortex core 
is responsible to the local Josephson coupling with a su- 
perconducting STM tip with the same symmetry. For 
the anti-parallel vortex state, as shown in Fig. 3(a), the 
STM experiments with an odd-frequency s-wave super- 
conducting tip may observe the spatial variation similar 
to the LDOS. As shown in Eq. however, the zero 
energy DOS does not contribute to even-frequency pair 
amplitudes, 7-±, triplet ( r j E = 0), which is not responsible 
to a STM with an even-frequency chiral p-wave super- 
conducting tip. This might potentially be a signature of 
the odd-frequency pairing. 

Before closing, we discuss candidates materials of chi- 
ral p-wave superconductor and supernuid. One of the 
experimentally accessible superconductor is Sr2Ru04-4S 
where surface Andrecv bound state specific to chiral 
p-wave symmetry has been observed^ Recently, there 



have been proposed several heterostructures topologi- 
cally equivalent to chiral p-wave superconductor, e.g., 
topological insulator/spin-singlet s-wave superconduc- 
tor heterostructures ^r— semiconductor/spin-singlet s- 
wave superconductor junctions with strong spin-orbit 
coupling The chiral superfluidity in Fermi gases 

near p-wave Feshbach resonances and spin-orbit coupled 
Fermi gases are also promising candidate systemsi 23 ' 49 i 50 
We hope our theoretical prediction will be experimentally 
verified in these systems. 



IV. SUMMARY 

In summary, we have examined the relation of the 
odd-frequency s-wave pair amplitude ^.triplet (r, E) and 
LDOS N(r, E) for the quantum limit of CdGM vortex 
bound states in chiral p-wave superconductors, based on 
the BdG equation. The relation of Eq. (flU|) holds in the 
Majorana zero modes of a single winding vortex. Fur- 
ther, we have confirmed that J 7 ,, triplet (r, E) and N(r,E) 
for the bound states have the same spatial structures even 
in finite E at the core of the vortex with w = — 1, and ap- 
proximately at chiral domain wall of the w = — 2 vortex. 
These spatial structures of the odd-frequency pair am- 
plitude give valuable information, when we identify the 
odd-frequency pair amplitude in future experiments and 
discuss the associated anomalous interference phenom- 
ena. 



ACKNOWLEDGMENTS 

We thank K. Machida for fruitful discussions. This 
work was supported by JSPS (No. 2074023303) and 
"Topological Quantum Phenomena" (No. 22103005) 
KAKENHI on innovation areas from MEXT. 



1 F.S. Bergeret, A.F. Volkov, and K.B. Efetov, Phys. Rev. 
Lett. 86, 4096 (2001). 

2 F.S. Bergeret, A.F. Volkov, and K.B. Efetov, Rev. Mod. 
Phys. 77, 1321 (2005). 

3 M. Eschrig, T. Lofwander, T. Champel, J. C. Cuevas, J. 
Kopu and G. Schon; J. Low Temp. Phys. 147, 457 (2007); 
M. Eschrig and T. Lofwander, Nature Phys. 4, 138 (2008). 

4 T. Yokoyama, Y. Tanaka, and A. A. Golubov, Phys. Rev. 
B 75, 094514 (2007). 

5 Y. Asano, Y. Tanaka, and A. A. Golubov, Phys. Rev. Lett. 
98, 107002 (2007); Y. Sawa, T. Yokoyama, Y. Tanaka, and 
A. A. Golubov, Phys. Rev. B 75, 134508 (2007). 

6 J. Linder, T. Yokoyama, A. Sudb0, and M. Eschrig, Phys. 
Rev. Lett. 102, 107008 (2009). 

7 Y. Tanaka and A. A. Golubov, Phys. Rev. Lett. 98, 037003 
(2007). 

8 Y. Tanaka, A. A. Golubov, S. Kashiwaya, and M. Ueda, 



Phys. Rev. Lett. 99, 037005 (2007); Y. Tanaka, Y. 
Tanuma, A. A. Golubov, Phys. Rev. B 76, 054522 (2007). 
9 Y. Asano, Y. Tanaka, A. A. Golubov, and S. Kashiwaya, 
Phys. Rev. Lett. 99, 067005 (2007). 

10 T. Yokoyama, Y. Tanaka and A. A. Golubov, Phys. Rev. 
B. 78, 012508 (2008). 

11 T. Yokoyama, M. Ichioka, and Y. Tanaka, J. Phys. Soc. 
Jpn. 79, 034702 (2010). 

12 Y. Tanaka, M. Sato, and N. Nagaosa, J. Phys. Soc. Jpn. 
81, 011013 (2012). 

13 Y. Tanaka and S. Kashiwaya, Phys. Rev. Lett. 74 3451 
(1995); S. Kashiwaya and Y. Tanaka, Rep. Prog. Phys. 
63 1641 (2000); Y. Asano, Y. Tanaka, and S. Kashiwaya: 
Phys. Rev. B 69 (2004) 134501 (2004). 

14 Y. Tanaka and S. Kashiwaya, Phys. Rev. B 70, 012507 
(2004); Y. Tanaka, S. Kashiwaya, and T. Yokoyama, ibid 
71, 094513 (2005); Y. Tanaka, Y. Asano, A. A. Golubov, 



6 



and S. Kashiwaya, Phys. Rev. B 72, 140503(R) (2005); 
Y. Asano, Y. Tanaka, and S. Kashiwaya, Phys. Rev. Lett. 
96, 097007 (2006); T. Yokoyama, Y. Tanaka, and N. Na- 
gaosa, Phys. Rev. Lett. 106, 246601 (2011); Y. Tanaka, Y. 
Nazarov, and S. Kashiwaya: Phys. Rev. Lett. 90 167003 
(2003). 

15 V.L. Berezinskii: JETP Lett. 20 (1974) 287. 

16 A. Balatsky and E. Abrahams, Phys. Rev. B 45, 13125 
(1992). 

17 C. Caroli, P.G. de Gennes, and J. Matricon, Phys. Lett. 9, 
307 (1964). 

18 N. Hayashi, T. Isoshima, M. Ichioka, and K. Machida, 
Phys. Rev. Lett. 80, 2921 (1998). 

19 S. Kaneko, M. Hafiz, E. Kakizaki, K. Matsuba, N, Nishida, 
T. Kawakami, T. Mizushima, and K. Machida, J. Phys. 
Soc. Jpn. 81, 063701 (2012). 

20 M. W. Zwierlein, J. R. Abo-Shaeer, A. Schirotzek, C. H. 
Schunck, and W. Ketterle, Nature (London) 435, 1047 
(2005). 

21 N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 

22 D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001). 

23 T. Mizushima, M. Ichioka, and K. Machida, Phys. Rev. 
Lett. 101, 150409 (2008). 

24 T. Mizushima and K. Machida, Phys. Rev. A 81, 053605 
(2010). 

25 S. Tewari, S. Das Sarma, and D.-H. Lee, Phys. Rev. Lett. 
99, 037001 (2007). 

26 V. Gurarie and L. Radzihovsky, Phys. Rev. B 75, 212509 
(2007). 

27 T. Kawakami, T. Mizushima, and K. Machida, J. Phys. 
Soc. Jpn. 80, 044603 (2011). 

28 S. Yasui, K. Itakura, and M. Nitta, Nucl. Phys. B 859, 
261 (2012). 

29 S. Higashitani, S. Matsuo, Y. Nagato, and K. Nagai, S. 
Murakawa, R.Nomura, and Y. Okuda, Phys. Rev. B 85, 
024524 (2012); S. Higashitani, Y. Nagato and K. Nagai, J. 
Low Temp. Phys. 155, 83 (2009). 

30 Y. Asano and Y. Tanaka, larXiv: 1204.42261 

31 M. Matsumoto and R. Heeb, Phys. Rev. B 65, 014504 
(2001). 

32 Y. Kato, J. Phys. Soc. Jpn. 69, 3378 (2000). 

33 Y. Kato and N. Hayashi, J. Phys. Soc. Jpn. 70, 3368 
(2001); ibid 71, 1721 (2002). 

34 y Tanuma, N. Hayashi, Y. Tanaka, and A. A. Golubov, 
Phys. Rev. Lett. 102, 117003 (2009). 



35 F. Gygi and M. Schluter, Phys. Rev. B 43, 7609 (1991). 

36 Y. Tanaka, H. Takayanagi, and A. Hasegawa, Solid State 
Commun. 85, 321 (1993). 

37 The symmetric relation for E-^—E holds in the imaginary 
part of the Green's function e~^ w+1 ' e J- 3 , while the real 
part is anti-symmetric. See for example, Y. Tanaka, Y. 
Asano, A. A. Golubov, and S. Kashiwaya, Phys. Rev. B 
72, 140503(R) (2005); Y. Tanaka and S. Kashiwaya, Phys. 
Rev. B 70, 012507 (2004) and Ref. B 

38 A. S. Mel'nikov and V. M. Vinokur, Nature (London) 415, 
60 (2002). 

39 JA. Sauls and M. Eschrig, New J. Phys. 11, 075008 (2009). 

40 A. P. Mackenzie and Y. Maeno, Rev. Mod. Phys. 75, 657 
(2003). 

41 S. Kashiwaya, H. Kashiwaya, H. Kambara, T. Furuta, H. 
Yaguchi, Y. Tanaka, and Y. Maeno, Phys. Rev. Lett. 107, 
077003 (2011). 

42 L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008). 

43 L. Fu and C. L. Kane, Phys. Rev. Lett. 102, 216403 (2009); 
A. R. Akhmerov, J. Nilsson, and C. W. J. Beenakker, Phys. 
Rev. Lett. 102, 216404 (2009); Y. Tanaka, T. Yokoyama, 
and N. Nagaosa, Phys. Rev. Lett. 103, 107002 (2009); J. 
Linder, Y. Tanaka, T. Yokoyama, A. Sudbo, and N. Na- 
gaosa, Phys. Rev. Lett. 104, 067001 (2010). 

44 K. T. Law, P. A. Lee, and T. K. Ng, Phys. Rev. Lett. 103, 
237001 (2009). 

45 M. Sato, Y. Takahashi, S. Fujimoto, Phys. Rev. Lett 103, 
020401 (2009); M. Sato, Y. Takahashi, and S. Fujimoto, 
Phys. Rev. B 82, 134521 (2010). 

46 J.D. Sau, R.M. Lutchyn, S. Tewari, and S. Das Sarma, 
Phys. Rev. Lett. 104, 040502 (2010); T. D. Stanescu, J. 
D. Sau, R. M. Lutchyn, and S. Das Sarma, Phys. Rev. B, 
81, 241310 (2010). 

47 J. Alicea, Phys. Rev. B 81, 125318 (2010). 

48 A. C. Potter and P. A. Lee, Phys. Rev. Lett., 105, 227003 
(2010); J. Linder and A. Sudbo, Phys. Rev. B 82, 085314 
(2010); A. Yamakage, Y. Tanaka and N. Nagaosa, Phys. 
Rev. Lett. 108, 087003 (2012). 

49 Y. Inada, M. Horikoshi, S. Nakajima, M. Kuwata- 
Gonokami, M. Ueda, and T. Mukaiyama, Phys. Rev. Lett. 
101, 100401 (2008) and references therein. 

50 P. Wang, Z.-Q. Yu, Z. Fu, J. Miao, L . Huang, S. Chai, H. 
Zhai, and J. Zhang, larXiv:1204.T887l 



